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In the presence of nonrelativistic supersymmetry, a sharp fermionic collective mode similar to the Goldstino 
mode in high-energy physics was proposed to be realized in Bose-Fermi mixtures. The Goldstino mode is 
relaxed (a.k.a. decays) if supersymmetry is explicitly broken, which can be revealed as the broadening of 
the corresponding spectral function. We find that the situation shares many similarities with the electron spin 
resonance in magnetic systems and adopt the well-known Kubo-Tomita theory to perform a general analysis of 
the spectral function lineshape broadening of the Goldstino mode. 


I. INTRODUCTION 

Supersymmetry (SUSY), the intriguing symmetry relating 
bosons and fermions, has been of strong interest in elementary 
particle physics after the first theoretical model in interacting 
quantum field theow with supersymmetry was constructed by 
Wess and Zumino 111]. It may play a fundamental role because 
the mathematical structure of supersymmetry called graded 
Lie superalgebra is the only one consistent with that in rel¬ 
ativistic quantum field theory 12]. The first realistic super- 
symmetric version of the standard model known as the min¬ 
imal supersymmetric standard model was also later proposed 
to solve the hierarchy problem about the large discrepancy be¬ 
tween aspects of the weak force and gravity a . Regardless of 
the intense theoretical works on supersymmetry in the particle 
physics, none of the superpartners, which have identical prop¬ 
erties except for opposite statistics, of any known elementary 
particle has been found in experiments so far. Even though 
the Large Hadron Collider (LHC) recently confirmed the dis¬ 
covery of the Higgs particle , there is still no evidence of 
any supersymmetric particle. 

Despite the difficulties in confirming the relativistic SUSY 
in high energy physics, there have been a few theoretical 
proposals on the spacetime SUSY in condensed matter and 
atomic systems, e.g. the possible emergence of spacetime 
SUSY at quantum critical points in lattice models 0,0]. simu¬ 
lation of the spacetime SUSY in optical lattices in atomic sys¬ 
tems 12, and the emergence of spacetime SUSY at the bound¬ 
ary of topological phases in condensed matter systems Il2 . 
On a different front, nonrelativistic SUSY (a Bose-Fermi 
symmetry unrelated to space-time symmetry) has attracted 
considerable interest due to the recent experimental progress 
in mixtures of ultra cold Bose and Fermi atoms UlUl^ . The 
Bose-Fermi mixtures provide an opportunity to realize and 
study supersymmetry. An ultra cold superstring model was 
constructed The physical behavior of an exactly 

soluble of one-dimensional Bose-Fermi mixtures was inves¬ 
tigated 123 - 1^ . A general formalism to study such super- 
symmetric systems based on coherent state path integral was 
also set up 12511 . 

In Ref. 1^ . one of the present authors studied general 
properties of supersymmetric Bose-Fermi mixtures, in which 
bosons and fermions are supersymmetric partners of each 


other. In the presence of time-reversal symmetry, supersym¬ 
metry is broken by a chemical potential difference between 
bosons and fermions, fif — fib = Afi > 0, which still keeps 
the canonical Hamiltonian, consisting of kinetic energy and 
the potential terms, supersymmetric. Such systems support 
a sharp fermionic collective excitation 122 similar to that in 
supersymmetric high-energy theories, referred to as the Gold¬ 
stino. The spectral function of the Goldstino mode has a sharp 
i5-function peak at zero momentum k = 0 and frequency 
u! = Afi. However, in experiments, it is a tremendous task 
to tune the parameters in such Bose-Fermi mixtures to keep 
the canonical Hamiltonian supersymmetric. The supersym¬ 
metry of the canonical Hamiltonian is expected to be broken 
and the sharp i5-function peak at a; = Afi of the spectral func¬ 
tion will be modified due to the relaxation of the Goldstino 
mode, which is the main focus in this work. 


Studying the relaxation of the Goldstino, a.k.a. the broad¬ 
ening of the sharp b-function peak at a; = Afi of the spec¬ 
tral function corresponding to the Goldstino mode, is not as 
straightforward as it seems, because the Goldstino mode itself 
is a highly nontrivial zero-momentum mode involving the lin¬ 
ear combination of the fermion-boson pair at all sites (or all 
momentum) in the system. The present problem we encounter 
mimics the electron spin resonance (ESR) with anisotropy 
in magnetic systems. Comparing the two problems, we are 
interested in the relaxation of the <7 = 0 Goldstino mode, 
Q ^ fj^j’ while the ESR measures the q = 0 mode, e.g. 
motion of total spin in Heisenberg anti-ferromagnets (AFM), 
S'+ = S^. We adopt the Kubo-Tomita (KT) theory Iz^ . 

the well-known theoretical approach for the ESR, to perform 
a general analysis of the relaxation of the Goldstino mode. 


The paper is organized as follows. In Sec.|II]we introduce 
the model system followed by performing the general analy¬ 
sis of the broadening effects due to supersymmetry breaking 
using KT theory. We conclude the work in Sec. [HI] We briefly 
review the KT theory in App. and provide a perturbative 
analysis in the non-interacting limit in Ann.lBlwhich is consis¬ 
tent with the results in the more general KT theory approach. 
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II. KUBO-TOMITA THEORY APPROACH TO THE 
RELAXATION OE GOLDSTINO MODE 


The model Hamiltonian is defined as H = Hq + H' = 
Hi + H 2 + H', where 


Hi — 

H2 = f + V = 

= -tHj2[b]b, + f]h]+uY: 

(jk) j 


H' = 5UY^ 




( 1 ) 




( 2 ) 

(3) 


The unperturbed Hamiltonian, Hq = Hi + H 2 , con¬ 
tains the chemical potential term, Hi, and the canonical 
Hamiltonian,iT 2 , at the SUSY point. H' represents the SUSY 
breaking term. In the Bose-Fermi mixtures, there are only two 
possible SUSY breaking terms-(i):the interaction strength dif¬ 
ference between the onsite boson-fermi interaction and the on¬ 
site boson-boson interaction and (ii):the difference between 
the hopping strengths of bosons and fermions. We consider 
the former case (i) first and the whole discussions below can 
also be applied to the latter case (ii), which will be discussed 
near the end of the paper. The Goldstino field is defined as. 


will show a i5-function peak located at w = B^- However, the 
presence of the anisotropy in the Heisenberg term H'^^^ will 
have effects on broadening the linewidth and shifting the peak 
center, 

Comparing the Hamiltonian in the Goldstino problem with 
that in the Heisenberg AFM, we can have almost a one-to-one 
correspondence between these two problems. The chemical 
potential term Hi in the Bose-Fermi mixtures is similar to 
the Zeeman energy H^ in the Heisenberg AFM, the canonical 
Hamiltonian H 2 consisting of kinetic and interacting terms is 
similar to the isotropic Heisenberg term Haf, and the super- 
symmetry breaking term H' is similar to the anisotropy term 
which significantly changes the spectral function line 

shape. 

Back to the Goldstino problem, we first notice that several 
commutators can be greatly simplified. 


H,,H2 


= 0 , 


H',Hi 


= 0 , 


( 10 ) 


Q,H2 



[Q, Hi] = (/i/ - fib)Q = A/iQ, 


( 11 ) 

( 12 ) 


while [H', H 2 ] and [Q, H'] remain complex and are not corn- 
mutable. 

In order to apply the KT theory, which is briefly reviewed 
in App.lAlwe introduce the interaction picture representation 
of Q as 




(4) 






(13) 


where N = Nf + Nf, is the total number of fermions and 
bosons. 

Before discussing how the supersymmetry breaking term 
affects the Goldstino spectral function lineshape, we first point 
out how the present model system mimics the ESR problem 
in the Heisenberg AFM with anisotropy. For a Heisenberg 
AFM in the Zeeman magnetic field with anisotropy in the 
Heisenberg term, which breaks the spin rotation symmetry, 
the Hamiltonian is Hafm = Haf + H^ + H' , with 


According to the KT theory, we can identity the Fourier com¬ 
ponent Q{uj = Ajj.) = Q. On the other hand, the interaction 
picture representation of H' is 

= = H'{uj = 0;t) = H'{t), (14) 

where we identity H'{uj = 0) = H'. The Goldstino field at 
time t, Q{t), can be expanded in powers of the perturbation as 
Q(t) = where the general expression of n-th 



(5) 

order of Q is 

j 

(6) 

« Jo Jo 1 


(7) 


{je) a=x,y,z 


where we introduce 


Q^^\t)-H'{ti),--- ,H'{tr,) 


ESR only probes the <7 = 0 total spin mode, 

and without anisotropy the equation of motion for the 

total spin can be solved exactly as 

dS+ 


dt 


= i [Haf + H,, 5 +] = -iB,S+ ( 8 ) 

^ 5+(f) = 5+(0)e"*^^‘. (9) 

Therefore, the spectral function corresponding to the retarded 
total spin Green’s function Gs{t) = —i6{t){[S'^{t), S'“(0)]) 


(15) 


[• • • [[A, B], C], • • • ,K]. The retarded Goldstino Green’s 
function can be expanded as 

00 

G«(i) = (0| {Q(f),Ql(0)} |0) = ^ G^(f),(16) 

n—0 

where 10) represents the (unperturbed) ground state. Accord¬ 
ing to Eqs. dnii-diiii, we can extract the Goldstino Green’s 
function up to any order and its corresponding spectral func¬ 
tion. 
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At zeroth order, we get 

G^{t) = = -i0(t)e-*^^‘.(17) 

Introducing Fourier transform, G^(a;) = f we 

get 


G^H = 


( 18 ) 


uj — Afi + i0+ 

and the spectral function at zero order is, = 

— (l/7r)/mG^(w), 


= (5(a; — A/r), (19) 


which shows a (5-peak located at A/r. The zeroth order (un¬ 
perturbed) result is consistent with that in Ref. 

For the first-order term, we have 

Gf (f) = dt' ({[Q, , Qt}) 

=f dt' {{[Q,H'],Q^]) 

* JO 

= -,e{t)e-'^'^\-it) {{[Q,H'],Q^}) , (20) 

where from first line to second line we use the fact that H 2 
commutes with Q. Using the identity {Q,Q^} = 1, we ob¬ 
tain that ({[Q, H'], Q1}) is real, which leads to the fact that 
the effect of the first-order term is only to shift the location 
of the (5-function peak of the spectral function, which will 
be presented below. For a more concrete result, we crudely 
mean-field decouple the expectation value, 

{{[Q,H'],Q^}) c -25C/(n5) ({Q,Qt}) = -2SUpU2l) 

where in the last equality we assume the density of the bosons 
is uniform and replace it by (rij) = N^/V = pb, with U being 
the whole system volume. After the mean-field decoupling, at 
first order we get the renormalized Goldstino Green’s function 

= G^{t) + Gf (f) ~ (1 -f i26Upbt) 

~ (22) 


where in the last line, we assume small SU and exponentiate 
the correcting term. In Fourier space, we get 


G«’(i)(a;) ~ 


Lo — (A/i — 2pb6U) + i0+ 
which gives the spectral function 

~ , 5 ( 0 ; - (Ap - 2pb5U)). 


(23) 


(24) 


The first-order perturbation indeed only shifts the (5-function 
peak. The result is consistent with the previous studies. 

In App. |Bj we perform the perturbative studies in the U = 0 
limit and we also obtain exactly the same result at first-order. 
For the second order term, the general form is 

= [ dti [\t2{{[Q\H'{G),H'{t2)],Q^)) 

Jo Jo 

= / dT{t-T){{[Q-H'{T),H'],Q^]). (25) 

^0 


If we follow KT theory to approximate 


= cr|/G(r), 

(26) 

al = {{[Q,H'],[H',Q^}), 

(27) 

({[g,7T'(T)],[iF',gt]}) 
({[g,iF'],[iT',gt]}) ’ 

(28) 


we can approximately rewrite G 2 if) as 

G^iO - -Goit)<Jo [ drit - T)/G(r). (29) 

Jo 

It is clear that (T^ cx (SU)'^ and /(r) is a “dimensionless” 
function which depends on variable r, and cr^ and /g(t) are 
both positive and real in this approximation. For the “time- 
independent” (Tq, we again crudely mean-field decouple it as 

al = ({[g, H '], [H', g^] }) iSU)^ {Qpl + 2pb) .(30) 

If we know the exact or approximate form of /g(t), we can 
get the exact or approximate form of G^(f). We do not find a 
way to get the exact form of /g(t), and below we take a plau¬ 
sible form of /g (t) jm, which was proposed in the studies in 
relaxation effects in nuclear magnetic resonance absorption, 
that leads to a reasonable result. 

Before discussing the general analysis of the second order 
term, we first point out that the main difficulty of extracting 
the exact form of /G('r) is that the time dependence of H'{t) 
can not be factored out similar to that for Qd^\T) due to the 
complicated form of [7T', H 2 \. Nevertheless, for a qualitative 
result in the “long-time” limit, t 00 , we can approximate 

nOO 

Gf (f) ~ -G^(t)cr|f / dT/cir) = -G^(<)cr^fT',(31) 
Jo 

where r' = dr fair). The renormalized Goldstino 
Green’s function containing up to second-order perturbation 
is 


G^d^)(l) _ Go(f) + Gi(f) + G2(f) 

~ G^(f) (1 -f i 2 pb 6 U - cFc^'t) 


which gives the spectral function 




TT [w — {Ap — 2pbSU)]^ + {(Tqt')'^ 


(33) 


We can see that the lineshape is a Lorentzian function with 
width Oqt' oc (SU)'^ centered around Ap — 2pbSU. (Jqt' 
corresponds to the “relaxation-rate” of the Goldstino or we 
can define a Goldstino “relaxation time” via the relation 
I/Tg = o'g''’^ O'- (dU)^. For further strengthening the con¬ 
jecture of the Lorentzian lineshape of the spectral function of 
Goldstino Green’s function, in App. |B} we provide the per¬ 
turbative calculations up to second-order of SU in the limit 
U = 0. We indeed find that the second-order self-energy 
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contains an imaginary part proportional to square of the per¬ 
turbation, oc (SU)'^, and the lineshape of the spec¬ 

tral function of Goldstino Green’s function is modified to be 
Lorentzian with width ~ oc {SUy, consistent with 

the results presented here using KT theory approach. 

In the general circumstance, it is difficult to get an analyt¬ 
ical result unless we know the exact form of /g(t)- Since 
/g(t) is a dimensionless function, the simplest form it can 
be, which is well-defined at any time, is 

/g(t) = (34) 


where tq is some “characteristic time” with unit of energy 
inverse (the same unit to time since we set H = 1). The 
simplest form of tq is a function consisting of the interac¬ 
tion strength inverse U~^, the hopping strength inverse 
and the chemical potentials of fermions and bosons, 

To = h{U ~^With the form of /g(t), we 
get 


G^{t) ~ -G^{t)(j}.Ty 
Therefore, 


e ^ 0 - 1 - 1 - 


(35) 




~ e‘^G^oG'^(i)e*2pi.-5c^e-‘"G^ot ^ 

n—0 


— ( 
nl V 


_ 2 2 -tTo\ 

^G^O 6 J ■ 


(36) 


The series converge very rapidly for ctqto < 1. The leading 
term with n = 0 gives the spectral function 




^2 ^2 
TT [W 


_ <^G^0 

(A/i - 2pbSU)f 



which is a Lorentzian function with half-width ctqTo ~ 

If acTo —>■ oo, which means tq —>■ oo since ctg Pb^U < 
1, Eq. drsl l can be Taylor expanded in powers of I/tq giving 
the leading term. 


G^{t)aGTo oi - G^{t)al:^—. 


(38) 


Note that in this limit, the leading second-order term is r- 
independent. If we ignore higher-order terms in the Taylor se¬ 
ries, the leading result corresponds to the approximation that 
— 0, which means the time-dependence of H'{t) 
is dropped. In this limit, the renormalized Goldstino Green’s 
function becomes 


G^[t)(TGTQ (39) 

which gives spectral function 

1 (Am-2p^5U)1=^ 

(uj)aGTo oi — . e . (40) 

\/2^ 


The spectral function lineshape becomes a “Gaussian distri¬ 
bution” with width (standard deviation) ctg cx \6U\, which 


is consistent with the first situation originally discussed by 
Kubo and Tomita and briefly reviewed in the App.lAl in which 
the Fourier components of perturbation, i.e. H'(uj), are time- 
independent. 

The crossover of the lineshape from a Lorentzian function 
to a Gaussian distribution can be possibly observed experi¬ 
mentally. In the present case, the commutator, [H' = 

\H', T] ^ 0, does not vanishes because the perturbation, H', 
is not commutable with the hopping term, T. In the cold 
atom systems, it is possible to gradually decrease the hopping 
strength by tuning the potential depth iGOl - l^ . which makes it 
possible to observe experimentally the change of the lineshape 
and the crossover from a Lorentzian with half-width oc (SU)^ 
to a Gaussian distribution with width oc \SU\ before entering 
the Mott-insulating phase. 

So far we have focused on the case in which the SUSY 
breaking term is the difference between the strength of the 
onsite fermion-boson and onsite boson-boson interaction. On 
the other hand, if we consider the other possible SUSY break¬ 
ing term, case (ii):the hopping strength difference between 
bosons and fermions, whose Hamiltonian is 

H'h=SthJ2fjfk, (41) 

0'fc> 

most of the discussions above can be directly applied to this 
case with H' H'^. For the first-order term, we find that 

,<9^}) = ((4) + (nO) = 

k 

where we introduce = — X]{e } with {e^} being 

the unit vectors that connect a site to its nearest neighbor sites. 
Following the previous discussions, we obtain, at the first or¬ 
der, the spectral function as ~ (5 (w — (A/r -f "f6th)), 

which shows the shift of the 5-peak. The second-order term 
involves 

({[0. , K.o'] }> = ((„*) + („/)) 

k 

= (43) 

where we introduce I]{e^} S{e,.} e*’' 

Following the previous discussions, we conclude that ctq ^ 
{SthY and the spectral function becomes a Lorentzian func¬ 
tion with half-width o-gTo ,th^ y y Pb^)^ 

where h'{x, y,z,w) isa phenomenology function that has the 
correct unit. Furthermore, if we experimentally “tune down” 
the interaction strength to make H 2 ] = [H/^y V] 0, the 
time-dependence of (f) can be dropped. The spectral func¬ 
tion lineshape in this limit becomes a Gaussian distribtution 
with width oc \Sth\, which can be experimentally confirmed. 


III. CONCLUSION 

We study the relaxation of Goldstino mode in Bose-Fermi 
mixtures due to the supersymmetry breaking. We adopt the 
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well-known Kubo-Tomita theory in electron spin resonance 
theory to find that the spectral function of Goldstino Green’s 
function, in general, is broadened to be a Lorentzian function 
with width proportional to the perturbation strength square. 
In the limit where either the hopping strength or the interac¬ 
tion strength vanishes, the lineshape becomes a Gaussian dis¬ 
tribution with width linearly proportional to the perturbation 
strength. 
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Appendix A: Brief Review of Kubo-Tomita Theory 


Most of the physical quantities can be obtained from the 
corresponding Green’s functions G{t), e.g. retarded Green’s 
function, real-time ordered Green’s function, imaginary-time 
Green’s function depending on the convention you use. The 
typical way to examine the resonance “frequency” is to ex¬ 
amine the corresponding Green’s function in Fourier space, 
G{u}) = J dujG{t)e^‘^*. In the noninteracting limit, H = Hq, 
the Fourier components give sharp resonance peaks at certain 
frequencies, G{uj) = G{uja)5{uj — uja), and this means 
the form of the Green’s function at time domain should be 
G{t) = In the presence of (weak) interaction 

H', the sharp resonance lines will get shifted and broadened. 
The conventional way to examine the effects of the interac¬ 
tion is to use equation of motion theory to expand the Green’s 
function in powers of t. Dehning the retarded Green’s func¬ 
tion G{t) = —i0(f)(O(f)O'l'(O))o, we can use equation of 
motion theory to expand 0{t) as 

t G 

0{t) = 0 ( 0 ) + -.[0{Q),H] + • ■ ■ + -[ 0 ( 0 ); H,--- ,H], 

% n\ 

(Al) 

where we introduce [A\ B,G, - ■ ■ ,K] = 

[■ ■ ■ [[A, B],G] ■ ■ ■ ,K]. The results above can be ob- 
tianed as follows. From Heisenberg equation of motion, we 
know 


dOjt) 

dt 




(A2) 


Integrating both sides with the initial condition 0{t = 0) = 
0(0), we obtain 


0{t) 


0 ( 0 ) + 



dti[0{ti),H]. 


(A3) 


By iteration, we get the result above. We then obtain 


G{t) = -ze{t) ^ ^ ^ I [0(0); if, • • • , iF], Ot (0)I ^ 

= Go{t) + Gi(t) + G2{t) + ■ ■ ■ + Gn{t)- (A4) 


Despite the simplicity of the expression above, it is, in most 
cases, difficult to analyze it analytically. The main reason is 
that we are interested in extracting the results perturbatively 
in powers of H', but, however, it is not easy to extract the 
perturbative results in powers of H' in this approach. Besides, 
if we go to the noninteracting limit, where the strength of the 
perturbation H' (x e vanishes, the perturbative result should 
go back to the noninteracting result as 

Hm,^oG{t) = Y, (A5) 

Ct 

In the presence of perturbation H', we may assume gener¬ 
ally that the coefficients of the exponential functions become 
function of t, 


G(e,f) = ^G„(e,f)e-*-“‘, (A6) 

06 


which gives the limits 

^oGa (c, f) — Go,, (A7) 

which is independent of t. 

For extracting the results perturbatively in powers of H', 
Kubo and Tomita (KT) introduced another way of performing 
the expansion in powers of the perturbation H'. Below, we 
briefly review the KT theory. 

As illustrated above, the most convenient way to obtain the 
expansion of G{t) is to solve the equation of motion, 

= [0{t), H] = [0(f), Ho + H’] . (A8) 

We are interested in performing the expansion in powers of 
if', and in order to achieve that, we switch to the interaction 
picture and introduce 

0{t) = (A9) 

where Hq is the unperturbed Hamiltonian. In the interaction 
picture, we obtain the equation of motion, 

= (AlO) 

where we also introduce ii'(f) = gi-Hot jj/^-iHot^ Integrating 
both sides of the equation with the initial condition O* {t = 
0) = 0(0) = O and we get 

0*{t) = 0+- f dh[0*{G),H'{-G)]. (All) 

* Jo 

By iteration, we get 

0*(f) = 0 +V(^)-" / / dtn-l-" dhx 

n^l -^0 -^0 

X ,ii'(-f„)](A12) 

Plugging the expansion above to Eq. ( IA9b gives 

0{t) = + ■ ■ ■ + + ■ ■ ■ (A13) 
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where the general expression of (t) is 







dti 




dtn X 




(A14) 


with 

O(0(i) = (A15) 

H\t) = (A16) 


The general way of evaluating the expression is to expand 
0(0) (f) and H'{t) in a Fourier series. We note that it is not 
always possible to decompose O^o) (t) and H'{t) into Fourier 
series exactly, but, fortunately, in most problems it is possible. 
There are two general situations depending on whether or not 
the Fourier components of and i.e., 0{uJa) 

and H'{uja), develop time dependence. Let us first focus 
on the first situation in which the Fourier components are 
“time-independent”, and the generalization to the second case 
can be straightforward carried over once we know all the 
logics developing in the first case. 


With the expansion of 0{t) in the powers of H' , the Green’s 
function can be straightforwardly expanded as 

G{t) = Go(f) + Gi(<) + G 2 {t) + • • • , (A17) 

where 

G„(f) = [ dti f dt2--- f dtn{-i)0{t) X 

Jo Jo Jo 

X ^|[0(°)(f); H'iG), • • • , Ot(0)|y 

(A18) 

The result we obtain above is from straightforward expansion 
in orders of the perturbation H' without making many 
assumptions, except for assuming H' can be treated as a 
perturbation. Hence, the general formula above should be 
valid in most (if not all) of the situations. 


In the first case, the 0^°) (t) and H'{t) can be fully decom¬ 
posed to a Fourier series. 


O0)(t)=^O(cc„)e-'-“‘, 

Ol 

(A19) 

Oi 

(A20) 

where the Fourier components O(uJa) and H'ioja) 
independent”. We can use the identity 

are “time 

e^Be — B + [A^ B\ + — [A, [A, B]] + — [A, [A, [A, H]]] 

+ ••• 

(A21) 

to get 


\0 (^U) ex') ^ Uq\ — CJq;0(cJq-), 

(A22) 

a 

a 

3 

II 

a 

(A23) 


Then we can repress Gn(t) as 


where 


Gn{t) 


a 


^ G '' ^ ^ (^1 ^1 ? ‘ ‘ ‘ ; (cCcK, , • • • 

Wl UJn 


• • • , — 


27ri 



e^*[p{p + iuJi){p + iuji + iuj2) ■ ■ ■ {p -f ^dp, 

e=i 


Fn{uia]UJl, ■ ■ ■ ,UJn) = ( ) [0{uJa)] ■ , H'(oJn)] , O') (0) 


(A24) 


(A25) 

(A26) 


We also know that O^t = 0) = and we can 

introduce 


,a;„) 



= ^F*{uJa:(^P\0Ji,--- ,U)n)- (A27) 

P 


From conservation of frequency, we know F* vanishes unless 


cUq. — u)p -f ^ ^ u)^ — 0. (A28) 

t=\ 

The reason that we introduce the function wi, • • • , uJn) 
is due to the fact the it is easier to analyze the function, 
whose result can be obtained by summing the residues of the 
simple poles located at p = 0, — iwi, — iuJi — iuj 2 , • • • , 
and etc. We can then write down the general expression along 
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with the constraint of Gn{t), 


metric terms is 


Gn{t) =”e 

a,0 Lj 

^^UJa—t^l3+J2e = l‘^t=0’ (A29) 

where we introduce the abbreviations, • • • J2uin’ 

gn { t -, Lo ) = g„(f; Wi, • • • ,Wn), and F*{uJa,uj/3;ui) = 

F*(w„, Wi, • • • ,u!n). Let’s try to examine some simple 
cases. For the first order perturbation, we have 


Gi(f) = -ie{t) ^ e— ^ -X 




a,p ci;i 

xF^{uJa,ujp;oji)Sui^ — U]I3+LJ1=0 

- itF^{uJa,u;a;0) + 


E l - g*(<^o-l^/3)t 

- F*{uJa,UJi3;-{u;a- UJfj)) 

^ Wq, — tU/J 


a . 

+E 


itF^ (cUq, , lUq^ , 0) -f 
{F^{uJa,U;p;-{uJa - W/j)) + 

, (A30) 


^13 

j ^/s)) 


where 


G 2 (f) cx—i0(t) ^2 e 


--F*{u;c.,ojc.;0,0)- 


-E 

7/0 


1 - iujjt - e 


.^2 (tOo, , 6 Uq. , , ^ 7 ) ■ 

(A34) 


Since we already know that the first-order perturbation only 
shift the location of the resonance line, let us see focus on the 
effects of second-order perturbation due to G 2 (f). Combining 
the zeroth order term with the second order term, we obtain 


G( 2 )(i)~Go(f) + G 2 (f) 


1 


i-E 


1 — iuj^t 


-cr^^ - -rcr. 


7/0 2 


aO 


-*^wEE 


Fpi^cn^a) ^-i(uia+uj-,)t ^2 


a 7/0 




a-y 


07 \ SF)^ - 


7/0 


-ie(t) J2 

a,77^0 


^07-^0 -i(ujc+uj-,)t 


w;: 


(A35) 


where we define = F 2 *(wa, Wq,; W 7 ,-w-),)/Fq (wq,, Wq,). 
The matrix element in the frequency space forms a trace and 
we can use the trace identity to rewrite 


F7a;o,a;o;0) = ^ | [0(w„), iT'(O)], (a;o)|^(A31) 

F*{u}a,t^l3]-{uJa “ UJ/j)) -f F* {oj fS] UJa] UJa “ U!p) 

= ^|[0(wo),i7'(w;3 -a;o)],Ot(u;;3)| -f 

+1 [O (t0;9), iT' (tOo - u;;3)], (wo) I y (A32) 

Using the fact H'"^{uj) = and similarly 0'^{uj) = 

0(—uj), we see that the terms in [...] are all real, which sim¬ 
ply indicates that the first order perturbation will only give 
resonance line “shift” in the Fourier space. The second-order 
perturbation contains more complex terms. The general ex¬ 
pression is 

G2{t) = ^g2{t-,w^,u;s) x 

oj 

xF^{uia,UJp\UJ^,UJs) X 6A33) 

There are several choices for satisfying the condition. 

(1) : UJa = w-y = 0 = (2) : Wa = W/3, w-y = 

^ 0. (3) . cUq. y u)^ — 0,w-y — (u^o ^z?)' 

(4) : Wo ^ Wy3, = 0,W5 = -(wo - ijp). (5) : = 

— {oJa — ^ 0. Below, we will only keep the terms 

with Wa = w^ and w-y = —W 5 and ignore other highly asym¬ 
metric terms with Wq, W/ 3 . The result ignoring highly asym- 


cr 


2 _ 
a7 — 


-F 2 *(Wa,Wa;W^,-W^) 

Fo*(Wa,Wa) 


(A36) 


[0(Wa),i7'(w-y)] , [i7'(-W.^),0'l'(Wa)] 


0 (Wa),Ot (Wo) 



Therefore, ij\^ and a^Q will be positive and real in this ap¬ 
proximation. The spectral function can be extracted according 
to A(w) = — (l/ 7 r)/mG(w), with G(w) = / dfG(f)e®“*, 

A( 2 )(w) =^F*(Wa,Wa) 




+EE Fp (Wa,Wa) ^ 

a 7/0 

Namely, the spectral function is the Gaussian line shifted by 
S 7/0 '^ 07/^7 from the original center w^ and broadened to 
width CTao- Besides, there is an additional resonance line, 
called a “satellite line” by Kubo and Tomita, with the rela¬ 
tive intensity at w^ -I- w-y. Therefore, we can see that 

if the perturbation can be completely decomposed into Fourier 


1 

\/27^ 
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series whose Fourier components do not develop time depen¬ 
dence, the perturbations broaden the original resonance line 
to a Gaussian line along with shifting the center, and there 
are additional resonance peaks (satellite lines). If the satellite 
lines are too close to the original resonance line peak, the per¬ 
turbation will break down. Mathematically the perturbation 
breaks down if lo^ <C <Ja-y or aao- 

Now let’s shift our focus on the second case in which the 
Fourier components of H'{u}a) develops time dependence. 
We will see that the Gaussian line will be modified to be 
a Lorentzian line, which was called motional narrowing by 
KT. In this case, the unperturbed Hamiltonian consists of two 
parts,iJo = Hi+ H 2 , which satisfy the conditions 


(A39) 

[i^2,O]=0. (A40) 

Hi is assumed to be not commutable with the operator O. 
With the conditions, we can straightforwardly adopt the ap¬ 
proach from the first case. In the interaction picture, the un¬ 
perturbed motion of the operator O can be expanded as 


q(0) _ 


Oi 

(A4I) 


where the Fourier component 0{u!a) satisfy 

[0(w„),iFi] =w„0(w«), (A42) 

[0(w„),iJ2] =0, (A43) 


Similarly, we also introduce 


H'it) 


^iHotjji_ ^iH2t^iHitjji 

= (A44) 

OL 


where we introduce 


^ (A45) 

a 

= H'{uJo^;t). (A46) 


With the new definitions, we can adopt all the approaches il¬ 
lustrated in the first case with the replacement of H'iuja) by 
the time-dependent H'{uia;t). The general term of the expan¬ 
sion is 

dti - ■■ f X 

Jo 

^ — - ^ 

X ^ I [©(wq,); iJ'(w-y; f 1), • • • , )], O'*' (w/j) I y 

(A47) 




Again, due to the conservation of the frequency , we require 


UJa - 


(A48) 


The zero-th order term is 

Go(f) = -z0(f)^e-*““y|oK),O+K)|y.A49) 
The first order term is 

Gi{t) f x 

a,13 7 

X (|| [C)(wc.), i?'i 1)], O'*(a;;3)I^ (A50) 

=-i6»(f) V X 

•’'0 

X ^ I [0(u;„), iF'(-; 0)], O* (w;3) I y (A51) 

which gives us the same result after integration to that in the 
first case, whose effects are to shift the center without broad¬ 
ening the resonance line. Let us focus on the second-order 
term, which is 

^ 2(0 = —iOit) X 

Q:,/3,7,(5 

^0 Jo 

X ^ I [O(Wa); iJ'(w^, f 1 - f2), '(W5,0)], O'* ) I y 

(A52) 

Again, we focus on the contributions from = oj^s, and 
iOj = —Lus and ignore highly asymmetric terms from uip 
Lop, which in most of the cases are not important. We get 

G2{t) ~ [ dh [' X 

a.7 “'0 “'0 

X ^ I [0(w„ ) ; ' (^7, f 1 - ^ 2 ), (-, 0)] , O'* (Wa ) I ^ 

= i0{t) V e-““‘ f drit - x 

«.7 •’'o 

X ^ I [ 0 (u;a); iT'(u; 7 , r), iT'(-W 7 , 0 )], O'* (Wa) I y 

(A53) 

Following Kubo and Tomita, we introduce 

^ I [0(c^„); H' (c^ 7 , r), iF' (-u; 7 ,0)], O* (cc„) | ^ 
^|o(cc„),OtK)|^ 

[0(t.„), if'K, r)], [i7'(-cc7,0), Ot^)] I ^ 
^|o(c^„),Ot(cc„)|^ 


U)^ + CJ-y + • ■ ■ + U)i/ - 0. 


(A54) 
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where has been introduced before and faj{T) is 
fa-fir) 

^ I [OK), il'K, r)], [H'i-io,, 0 ), Ot (cc„)] I ^ 

^I [OK), ii'K. 0 )], [H'i-uj^, 0 ), otK)] 

(A55) 

which is reduced to be 1 if becomes time indepen¬ 

dent. Then the expression for the second ordre term is 


G2{t) ~z0(f)^e-*‘"“‘^|oK),Onwc)|^ X 

X [ dT{t - (A56) 

Jo 

Now, let us combine the terms up to second order. Since first 
order term is the same to that in the first case, it will simply 
shift the resonance line, so we ignore it. Up to second order, 
we get 

G( 2 )(f)~Go(f) + G 2 (f) 

- G'o(f) X! (l “ ~ 


- Goit) 


K p~fo 


7 


Focusing on the simplest case with 7 = 0, we have 
~ Go(f)e"‘"“ 0 '^o 

At the long time limit, t 00 , we can approximate 
G( 2 )(f ^ 00 ) - Goit)e-'^"°ff°drtfMr) 

= Go(f)e-‘"“«‘("'+*"''), 

where we introduce r' = Re{fao{T)}dT and t" = 
/m{/ao('r)}(iT. Then the spectral function becomes 

A( 2 )K:x 1 - 


(A57) 

(A58) 

(A59) 

(A60) 


TT [w - K - <rMr")Y + 
where we can see that the line shape becomes Lorentzian! 
Note that if we set faoir) = fao^ independent of time, the 
shape will go back to Gaussian. In general, if we can get 
the “exact” form of fa-fir), we can obtain the exact form of 
the resonance line, which is, unfortunate, not possible in most 
cases. 


Appendix B: Perturbative studies in the 17 = 0 limit 


We present the perturbative studies of the broadening of the 
spectral function of Goldstino mode at G = 0. The Hamilto¬ 
nian is 77 = Hq + H' with 


Ho = -th ^ + fjfk 

O'fc) 


M/A ~ fJ-bNb, (Bl) 


H' = 

3 


(B2) 



■> 

(a) 



FIG. 1. (Color Online) (a) The Goldstino Green's function con¬ 
sisting of a fermion Green’s function (red dashed line) and a bo¬ 
son Green’s function (blue solid line) whose momenta are explic¬ 
itly summed, (b) The two-fermion-two-boson vertex arising from 

H = UbF 'djfl'j = X)p,q,k /p/p+q-k&q^k- 


P 

- 

(a) 



> = > 


(c) 


FIG. 2. (Color Online) First-order Feynman diagrams contributing to 
renormalize the Goldstino Green’s function, (a) A boson bubble in 
the fermion Green’s function sector renormalizing fif. (b) A fermion 
buble in the boson Green’s function sector renormalizing fib- (c) A 
two-fermi-two-bose vertex contributing to the first-order self-energy. 


For convenience, we introduce the imaginary-time-ordered 
Goldstino Green’s function 


a(r) = -(T. [Q(r)gt(o)])^ 

where we remind that the Goldstino field is 

The general expression for the perturbation is 

(-ir 


=- 1 ] 

m —0 

X A, 


m! 


d-Ti ■ ■ 


dTm X 


[H'(ri)..K'K)g(T)Qt(o) 


(B3) 


(B4) 


OfConn 


05) 


where the subscript 0 means the unperturbed ground state and 
conn means the connected Feynman diagrams. 

At zeroth order, we know that 

C;W(r)=-(T.[g(r)gt(0)]) 

= [fli'^)Mr)fpi0)bli0)]) 

p,q 

= [fpiO)flir)]){Tr [6,(r)6j(0)]) 

p.q 

= P)Gbir, p), (B 6 ) 

















gives (in imaginary time domain) 
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> ! > | > 

> '">'■> 

(a) (b) 


->- 

(c) 


FIG. 3. (Color Online) The second-order irreducible Feynman dia¬ 
grams. (a) Two cross interaction vertices, (b) A fermionic particle- 
hole bubble in the boson Green’s function sector, (c) A bosonic 
particle-hole bubble in the fermion Green’s function sector. 


where p) are the imaginary-time ordered 

fermion/boson Green’s functions, and we know 


Si{t) = I dTi X 


N 


Pb 


^Gfiri -T,p)gf{-Ti,p)gb{T,p) + (Bii) 


+PF'^gf{-T,p)gb{Ti,p)gb{T - ti , p ) - 


{ T -- Tl , q ) gf {- Tl , p)gb (ti , p) . 
p,q -* 

(B12) 


Going to Fourier space, we get 


p)= [6>(r)n/(ej^) - 6>(-r) (1 - n/(ej^))] e"p^, (B7) gi{iu;n) = g^°\iujn) {-PbUbf + PfUbf - pUbf) g^°\iuJn) 
gb{T, p) = — [9{t) (l -F nf,(ep)) + 0(— r)nt,(ey (B8) = {—2pbUbf) (B13) 


where we introduce = (e^^±l)“^ andcp^^ = ujp^^— 

p f /b, with being the dispersions of hopping. Introducing 

the Fourier transform C/(iw„) = drC/(T)e*“"'’’ with = 

(2n + 1)^f [3{2mrI(3) for fermions/bosons, we get 





p 


»/(ep) + ^fc(£p) 

iOJn + ep - Ep 


1 

iuin — Ap 


,(B9) 


where we note that ojn = {2n + \)f/(3 since the Goldstino 
field is a fermion. The zeroth order (unperturbed) retarded 
Goldstino Green’s function can be obtained as 

^ w + ^0+) =-(BIO) 

which is the same to that derived from the equation of motion 
presented in the main texts. For clarity in presenting the per¬ 
turbative studies below, we introduce the diagram presentation 
for the Goldstino Green’s function in Fig.flja), where the up¬ 
per dashed red line represents the fermion line and the bottom 
blue solid line represents the boson line and the momenta of 
fermion and boson lines are explicitly summed. Fig.lTfb) rep¬ 
resents the two-fermion-two-boson vertex. The arrows repre¬ 
sent the directions of the particles. 

For the first order term, we find there are three possible 
contributions as shown in Fig. |2] The first-order perturbation 


Therefore, we can see if we only include the first order term, 
the renormalized imaginary-time Goldstino Green’s function 
at first order is 


g^^\iuJn) - TIT -^- J-j -T, 

- (A/r - 2pbUbf) 


(B14) 


which gives 




1 

w — (Ap — 2pbUbf) + *0+ ’ 


(B15) 


which again is consistent with the result obtained by Kubo- 
Tomita theory presented in the main texts and leads to the 
spectral function 

A^^\u:) = S{u;-{Ap-2pbUbf))- (B16) 


For the second order, we find there are three contributions 
illustrated schematically in Fig. [3] The derivation for the 
second-order terms are straightforward (but tedious), and we 
present the final results below. Fig [3^a)-(c) contributing to 
renormalize the imaginary-time Goldstino Green’s functions 

can be concisely written as g^°\ 

where the subscript (2) means the second-order and the sub¬ 
scripts 1 — 3 label the contributions from Fig. [3ja)-(c). The 
self-energies E contributed from each term are 

















11 


(2)r. . Ulp V- ( + ^biei)] ^^/(ep+q-k)-»^/(ep-4) rifie{) + nbi-el 

^1 2^ I + g/ + g/ _ g& _ / 

p,q,k ^ i- tp + eq ej^ S+q-k 


£p + 4 + ^p+q-k) 


K(4) - «b(ep)] 

iuJn + Cp + eq - - 

”/(ep+q-k) +”6(4 - 4) 

p+q-k 

^^6(4) + ”/(-4 + 4 - 4+q-k) 

K(ep - ’^b(4)] 

iuJn -e^-e6-Fe6- 

”/(4+q-k)+”6(4"4) 

4+q-k 

”/(4) “ ”/(-4 + 4 + 4+q-k) 

h/(4) + ”6(4)] 

H- 

S ^ 

”/(4+q-k) -”/(4 -4) 

>+q-k 

nh(e^) - nb{-e{, + ei + 4+q-k) 


2 '' _/Vl /2 


E 

p,q.k 


:b _ fb _ cl 

q ^k ^p+q-k 


-n/(e/) n/(e^+k-p)-«/(4-ep) ’^/K) ++ 4 + ^q+k-p) 


iwn - - 4 + ep + Cq+k-p 


+ - 


»^b(ep) »^/(eq+k-p)+^b(4 - 4) ^/(eq)-^/(-4 + ep + 


^q+k-p-^ 


^q + ^k ^p ^q+k-p 


:3i7) 


(B18) 




p,q.k 


nh(e^) rb/(e{) + rbh(ep + eq+k-p) ’^b(ek) + ’^/( 


i-.{ 


""q+k 


-p) 


fb 

q+k-p 


/(ej^) ’^/(e^)-»^/(e^ + eq+k-p) ^f>(ek) “’^b(-eq + ^p + ^q+k-p) 


- fb — cl _ fb 

''k '^P '"q+k-p 


(B19) 


The second-order self-energies in the Matsubara frequency 
domain are very complicated. However, since we are only 
interested in the imaginary part of the self-energy in the real 
frequency domain obtained by iujn —?► w -f iO’*' and extract its 
qualitative behavior. In this way, we can see that 2 3 ('^) 
all contains imaginary parts which are proportional to square 
of the perturbation strength, i. e. {uj) = oc Ugp. 

In the presence of the imaginary part of the self-energy, the 
retarded Goldstino Green’s function becomes (we ignore the 
real part of since it only contributes to shift the loca¬ 

tion of the peak) 




1 

w — (Ap — 2pbUbf) — ir( 2 ) 


(B20) 


which leads to the spectral function 




r( 2 ) 

(w-(Ap-2pBC/BF))2 + (r(2))2’ 


which is a Lorentzian function with width ~ F^^^ oc Upp- 
Therefore, we can see that the perturbative results at U = 0 
limit are qualitatively consistent with the results obtained by 
Kubo-Tomita theory. 
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